%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plots for skillbias_main
% this version: October 2010
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function plots_pub(Respmed,Resplow,Resphigh,spec,HORIZON_impresp,nvars,sample,rw,prem_choice,prod_choice,price_choice,hours_choice,do_break,break_choice,...
    prior,nlagsvar,decay,pndx,hndx,spndx,pricendx,wLndx,relendx,wHndx,hHndx,hLndx,do_level,do_sign)

tt = (0:HORIZON_impresp);
zerobar = zeros(1,HORIZON_impresp);

% Gali plots
if spec == 1 || spec == 2 || spec == 3
    
    figure(1);
    set(gcf, 'PaperSize', [13 3]);
    set(gcf, 'PaperUnits', 'centimeters');
    set(gcf, 'PaperPositionMode', 'manual');
    set(gcf, 'PaperPosition', [0 0  13 3]);
    subplot(1,4,1); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(pndx,:)],'-k',tt,[nan,Resplow(pndx,:)],'--r',tt,[nan,Resphigh(pndx,:)],'--r')
    if prod_choice == 1
        %axis([0 10 -0.1 3.9]);
        set(gca,'YTick',[0 1.5 3])
        ylabel('TFP measure');
    else
        axis([0 10 -0.1 1.0]);
        if spec == 1
            if sample == 6
                if prior == 1
                    ylabel('postwar bivariate');
                elseif prior == 2
                    ylabel('Gali');
                    title('productivity');
                end
            elseif sample == 7
                ylabel('1955-2000')
            else
                ylabel('recent bivariate');
            end
        elseif spec == 2
            ylabel('baseline');
        else
            ylabel('baseline');
        end
    end
    %     xlabel('Quarters');
    set(gca,'FontSize',7);
    subplot(1,4,2); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(hndx,:)],'-k',tt,[nan,Resplow(hndx,:)],'--r',tt,[nan,Resphigh(hndx,:)],'--r')
    if prod_choice == 0 && spec == 1 && sample == 6 && prior == 2
        title('tot. hours')
    end
    %     xlabel('Quarters');
    axis([0 10 -1 0.5]);
    set(gca,'FontSize',7);
    if spec == 1 || spec == 23
        subplot(1,4,3); plot(tt,[0,zerobar],'-k',tt,[nan,zeros(1,size(Respmed,2))],'-k')
    else
        subplot(1,4,3); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(spndx,:)],'-k',tt,[nan,Resplow(spndx,:)],'--r',tt,[nan,Resphigh(spndx,:)],'--r')
    end
    if prod_choice == 0 && spec == 1 && sample == 6 && prior == 2
        title('premium')
    end
    axis([0 10 -0.2 0.2]);
    %     xlabel('Quarters');
    set(gca,'FontSize',7);
    if spec == 1 || spec == 23
        subplot(1,4,4); plot(tt,[0,zerobar],'-k',tt,[nan,zeros(1,size(Respmed,2))],'-k')
    else
        subplot(1,4,4); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(relendx,:)],'-k',tt,[nan,Resplow(relendx,:)],'--r',tt,[nan,Resphigh(relendx,:)],'--r')
    end
    if prod_choice == 0 && spec == 1 && sample == 6 && prior == 2
        title('rel. hours')
    end
    set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 3.5])
    %     xlabel('Quarters');
    axis([0 10 -0.7 0]);
    set(gca,'FontSize',7);
    if prod_choice == 1
        print('-depsc2',['gali_BFK_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
    else
        if spec == 1
            print('-depsc2',['gali_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
        elseif spec == 2
            print('-depsc2',['galipluspremium_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
        elseif spec == 3
            if sample == 3 || sample == 4
                print('-depsc2',['galiplusprice_spec_' int2str(spec) '_sample' int2str(sample) '_rw' int2str(rw) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
            else
                print('-depsc2',['galiplusprice_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
            end
        end
    end
    
    % Fisher plots
elseif spec == 5
    
    figure(1);
    subplot(2,4,1); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(pricendx,:)],'-k',tt,[nan,-Resplow(pricendx,:)],'--r',tt,[nan,-Resphigh(pricendx,:)],'--r')
    title('price');
    ylabel('IBT shock');
    axis([0 10 -1.5 1]);
    subplot(2,4,2); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(pndx,:)],'-k',tt,[nan,Resplow(pndx,:)],'--r',tt,[nan,Resphigh(pndx,:)],'--r')
    title('prod.')
    axis([0 10 -0.1 0.6]);
    subplot(2,4,3); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(hndx,:)],'-k',tt,[nan,Resplow(hndx,:)],'--r',tt,[nan,Resphigh(hndx,:)],'--r')
    title('tot. hours');
    axis([0 10 -1 1]);
    subplot(2,4,4); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(spndx,:)],'-k',tt,[nan,Resplow(spndx,:)],'--r',tt,[nan,Resphigh(spndx,:)],'--r')
    title('premium');
    axis([0 10 -0.2 0.3]);
    subplot(2,4,5); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(nvars+pricendx,:)],'-k',tt,[nan,-Resplow(nvars+pricendx,:)],'--r',tt,[nan,-Resphigh(nvars+pricendx,:)],'--r')
    ylabel('INT shock');
    axis([0 10 -1.5 1]);
    xlabel('Quarters');
    subplot(2,4,6); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pndx,:)],'-k',tt,[nan,Resplow(nvars+pndx,:)],'--r',tt,[nan,Resphigh(nvars+pndx,:)],'--r')
    axis([0 10 -0.1 0.6]);
    xlabel('Quarters');
    subplot(2,4,7); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hndx,:)],'-k',tt,[nan,Resplow(nvars+hndx,:)],'--r',tt,[nan,Resphigh(nvars+hndx,:)],'--r')
    axis([0 10 -1 1]);
    xlabel('Quarters');
    subplot(2,4,8); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
    axis([0 10 -0.2 0.3]);
    xlabel('Quarters');
    set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 8])
    if sample == 3 || sample ==4
        print('-depsc2',['fisher_premium_spec_' int2str(spec) '_sample' int2str(sample) '_rw' int2str(rw) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
    else
        print('-depsc2',['fisher_premium_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
    end
    
    % SBT plots
elseif spec == 6 || spec == 9 || spec >= 25 && spec <= 31
    if do_sign == 1
        
        if spec == 9
            figure(1);
            subplot(1,4,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(spndx,:)],'-k',tt,[nan,Resplow(spndx,:)],'--r',tt,[nan,Resphigh(spndx,:)],'--r')
            axis([1 10 -0.5 0.5]);
            title('premium');
            ylabel('SBT shock');
            subplot(1,4,2)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(relendx,:)],'-k',tt,[nan,Resplow(relendx,:)],'--r',tt,[nan,Resphigh(relendx,:)],'--r')
            axis([1 10 -1.2 1.2]);
            title('rel. hours')
            subplot(1,4,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(pndx,:)],'-k',tt,[nan,Resplow(pndx,:)],'--r',tt,[nan,Resphigh(pndx,:)],'--r')
            axis([1 10 -0.5 0.6]);
            title('prod.');
            subplot(1,4,4)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(hndx,:)],'-k',tt,[nan,Resplow(hndx,:)],'--r',tt,[nan,Resphigh(hndx,:)],'--r')
            axis([1 10 -1 0.5]);
            title('tot. hours');
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 3.5])
            if sample == 3 || sample ==4
                display('no plot for this sample')
            else
                print('-depsc2',['SBT_SBTshock_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            end
            figure(2);
            subplot(1,4,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
            %                 title('premium');
            axis([1 10 -0.5 0.5]);
            ylabel('other T shock');
            subplot(1,4,2)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+relendx,:)],'-k',tt,[nan,Resplow(nvars+relendx,:)],'--r',tt,[nan,Resphigh(nvars+relendx,:)],'--r')
            %                 title('rel. emp.')
            axis([1 10 -1.2 1.2]);
            subplot(1,4,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pndx,:)],'-k',tt,[nan,Resplow(nvars+pndx,:)],'--r',tt,[nan,Resphigh(nvars+pndx,:)],'--r')
            %                 title('productivity');
            axis([1 10 -0.5 0.6]);
            subplot(1,4,4)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hndx,:)],'-k',tt,[nan,Resplow(nvars+hndx,:)],'--r',tt,[nan,Resphigh(nvars+hndx,:)],'--r')
            %                 title('hours');
            axis([1 10 -1 0.5]);
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 3])
            if sample == 3 || sample ==4
                display('no plot for this sample')
            else
                print('-depsc2',['SBT_otherTshock_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            end
        end
        
    elseif do_sign == 2
        if spec == 29 || spec == 30 || spec == 31
            
            figure(1);
            subplot(3,4,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(spndx,:)],'-k',tt,[nan,-Resplow(spndx,:)],'--r',tt,[nan,-Resphigh(spndx,:)],'--r')
            axis([1 10 -0.5 0.5]);
            title('premium');
            ylabel('supply shock');
            %axis([1 10 -0.4 0.6]);
            subplot(3,4,2)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(pndx,:)],'-k',tt,[nan,-Resplow(pndx,:)],'--r',tt,[nan,-Resphigh(pndx,:)],'--r')
            title('productivity');
            %axis([1 10 -0.1 0.6]);
            subplot(3,4,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(hndx,:)],'-k',tt,[nan,-Resplow(hndx,:)],'--r',tt,[nan,-Resphigh(hndx,:)],'--r')
            title('hours');
            %axis([1 10 -1 0.5]);
            subplot(3,4,4); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(relendx,:)],'-k',tt,[nan,-Resplow(relendx,:)],'--r',tt,[nan,-Resphigh(relendx,:)],'--r')
            title('rel. emp.')
            subplot(3,4,5)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
            ylabel('SBT shock');
            %axis([1 10 -0.4 0.6]);
            subplot(3,4,6)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pndx,:)],'-k',tt,[nan,Resplow(nvars+pndx,:)],'--r',tt,[nan,Resphigh(nvars+pndx,:)],'--r')
            %axis([1 10 -0.1 0.6]);
            subplot(3,4,7)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hndx,:)],'-k',tt,[nan,Resplow(nvars+hndx,:)],'--r',tt,[nan,Resphigh(nvars+hndx,:)],'--r')
            %axis([1 10 -1 0.5]);
            subplot(3,4,8); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+relendx,:)],'-k',tt,[nan,Resplow(nvars+relendx,:)],'--r',tt,[nan,Resphigh(nvars+relendx,:)],'--r')
            subplot(3,4,9)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+spndx,:)],'-k',tt,[nan,-Resplow(2*nvars+spndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+spndx,:)],'--r')
            ylabel('other T shock');
            xlabel('Quarters');
            %axis([1 10 -0.4 0.6]);
            subplot(3,4,10)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+pndx,:)],'-k',tt,[nan,-Resplow(2*nvars+pndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+pndx,:)],'--r')
            xlabel('Quarters');
            subplot(3,4,11)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+hndx,:)],'-k',tt,[nan,-Resplow(2*nvars+hndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+hndx,:)],'--r')
            xlabel('Quarters');
            %axis([1 10 -1 0.5]);
            subplot(3,4,12); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+relendx,:)],'-k',tt,[nan,-Resplow(2*nvars+relendx,:)],'--r',tt,[nan,-Resphigh(2*nvars+relendx,:)],'--r')
            xlabel('Quarters')
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 8])
            if sample == 3 || sample == 4
                display('not saved: choose other sample');
            else
                print('-depsc2',['SBT_tri_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            end
            
        elseif spec == 25
            
            figure(1);
            subplot(1,4,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(spndx,:)],'-k',tt,[nan,-Resplow(spndx,:)],'--r',tt,[nan,-Resphigh(spndx,:)],'--r')
            axis([1 10 -0.5 0.5]);
            ylabel('supply shock');
            xlabel('Quarters');
            subplot(1,4,2)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(relendx,:)],'-k',tt,[nan,-Resplow(relendx,:)],'--r',tt,[nan,-Resphigh(relendx,:)],'--r')
            axis([1 10 -1.2 1.2]);
            xlabel('Quarters');
            subplot(1,4,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(pndx,:)],'-k',tt,[nan,-Resplow(pndx,:)],'--r',tt,[nan,-Resphigh(pndx,:)],'--r')
            axis([1 10 -0.5 0.6]);
            xlabel('Quarters');
            subplot(1,4,4);
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(hndx,:)],'-k',tt,[nan,-Resplow(hndx,:)],'--r',tt,[nan,-Resphigh(hndx,:)],'--r')
            axis([1 10 -1 0.5]);
            xlabel('Quarters');
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 3.5])
            if sample == 3 || sample == 4
                display('not saved: choose other sample');
            else
                print('-depsc2',['SBT_trisupshock_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            end
            
            figure(2);
            subplot(1,4,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
            axis([1 10 -0.5 0.5]);
            ylabel('SBT shock');
            subplot(1,4,2)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+relendx,:)],'-k',tt,[nan,Resplow(nvars+relendx,:)],'--r',tt,[nan,Resphigh(nvars+relendx,:)],'--r')
            axis([1 10 -1.2 1.2]);
            subplot(1,4,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pndx,:)],'-k',tt,[nan,Resplow(nvars+pndx,:)],'--r',tt,[nan,Resphigh(nvars+pndx,:)],'--r')
            axis([1 10 -0.5 0.6]);
            subplot(1,4,4);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hndx,:)],'-k',tt,[nan,Resplow(nvars+hndx,:)],'--r',tt,[nan,Resphigh(nvars+hndx,:)],'--r')
            axis([1 10 -1 0.5]);
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 3])
            if sample == 3 || sample == 4
                display('not saved: choose other sample');
            else
                print('-depsc2',['SBT_triSBTshock_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            end
            
            figure(3);
            subplot(1,4,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+spndx,:)],'-k',tt,[nan,-Resplow(2*nvars+spndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+spndx,:)],'--r')
            axis([1 10 -0.5 0.5]);
            ylabel('other T shock');
            %             xlabel('Quarters');
            subplot(1,4,2)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+relendx,:)],'-k',tt,[nan,-Resplow(2*nvars+relendx,:)],'--r',tt,[nan,-Resphigh(2*nvars+relendx,:)],'--r')
            axis([1 10 -1.2 1.2]);
            %             xlabel('Quarters');
            subplot(1,4,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+pndx,:)],'-k',tt,[nan,-Resplow(2*nvars+pndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+pndx,:)],'--r')
            axis([1 10 -0.5 0.6]);
            %             xlabel('Quarters');
            subplot(1,4,4);
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+hndx,:)],'-k',tt,[nan,-Resplow(2*nvars+hndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+hndx,:)],'--r')
            axis([1 10 -1 0.5]);
            %             xlabel('Quarters');
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 3])
            if sample == 3 || sample == 4
                display('not saved: choose other sample');
            else
                print('-depsc2',['SBT_trioTshock_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            end
            
            figure(4);
            subplot(1,3,1); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+relendx,:)],'-k',tt,[nan,Resplow(nvars+relendx,:)],'--r',tt,[nan,Resphigh(nvars+relendx,:)],'--r')
            axis([1 10 -0.1 1]);
            title('rel. hours');
            xlabel('Quarters');
            subplot(1,3,2); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hHndx,:)],'-k',tt,[nan,Resplow(nvars+hHndx,:)],'--r',tt,[nan,Resphigh(nvars+hHndx,:)],'--r')
            axis([1 10 -1.5 0.5]);
            title('hours high');
            xlabel('Quarters');
            subplot(1,3,3); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hLndx,:)],'-k',tt,[nan,Resplow(nvars+hLndx,:)],'--r',tt,[nan,Resphigh(nvars+hLndx,:)],'--r')
            axis([1 10 -1.5 0.5]);
            title('hours low');
            xlabel('Quarters');
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 4.5])
            print('-depsc2',['SBT_triSBTshock_hHhL_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
            
            
            
        elseif spec == 26
            
            
            figure(1);
            subplot(1,3,1); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
            title('premium');
            axis([1 10 -0.1 0.6]);
            subplot(1,3,2); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+wHndx,:)],'-k',tt,[nan,Resplow(nvars+wHndx,:)],'--r',tt,[nan,Resphigh(nvars+wHndx,:)],'--r')
            title('wage high');
            axis([1 10 -0.1 0.6]);
            subplot(1,3,3); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+wLndx,:)],'-k',tt,[nan,Resplow(nvars+wLndx,:)],'--r',tt,[nan,Resphigh(nvars+wLndx,:)],'--r')
            title('wage low');
            axis([1 10 -0.1 0.6]);
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 4])
            print('-depsc2',['SBT_triSBTshock_wHwL_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '.eps'])
            
        elseif spec == 27
            
            relsndx = 6;
            figure(1);
            subplot(3,5,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(spndx,:)],'-k',tt,[nan,-Resplow(spndx,:)],'--r',tt,[nan,-Resphigh(spndx,:)],'--r')
            axis([1 10 -0.5 0.5]);
            title('premium');
            ylabel('supply shock');
            %axis([1 10 -0.4 0.6]);
            subplot(3,5,2)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(pndx,:)],'-k',tt,[nan,-Resplow(pndx,:)],'--r',tt,[nan,-Resphigh(pndx,:)],'--r')
            title('productivity');
            %axis([1 10 -0.1 0.6]);
            subplot(3,5,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(hndx,:)],'-k',tt,[nan,-Resplow(hndx,:)],'--r',tt,[nan,-Resphigh(hndx,:)],'--r')
            title('hours');
            %axis([1 10 -1 0.5]);
            subplot(3,5,4); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(relendx,:)],'-k',tt,[nan,-Resplow(relendx,:)],'--r',tt,[nan,-Resphigh(relendx,:)],'--r')
            title('rel. emp.')
            subplot(3,5,5); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(relsndx,:)],'-k',tt,[nan,-Resplow(relsndx,:)],'--r',tt,[nan,-Resphigh(relsndx,:)],'--r')
            title('rel. supply')
            subplot(3,5,6)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
            ylabel('SBT shock');
            %axis([1 10 -0.4 0.6]);
            subplot(3,5,7)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pndx,:)],'-k',tt,[nan,Resplow(nvars+pndx,:)],'--r',tt,[nan,Resphigh(nvars+pndx,:)],'--r')
            %axis([1 10 -0.1 0.6]);
            subplot(3,5,8)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hndx,:)],'-k',tt,[nan,Resplow(nvars+hndx,:)],'--r',tt,[nan,Resphigh(nvars+hndx,:)],'--r')
            %axis([1 10 -1 0.5]);
            subplot(3,5,9);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+relendx,:)],'-k',tt,[nan,Resplow(nvars+relendx,:)],'--r',tt,[nan,Resphigh(nvars+relendx,:)],'--r')
            subplot(3,5,10);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+relsndx,:)],'-k',tt,[nan,Resplow(nvars+relsndx,:)],'--r',tt,[nan,Resphigh(nvars+relsndx,:)],'--r')
            subplot(3,5,11)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+spndx,:)],'-k',tt,[nan,-Resplow(2*nvars+spndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+spndx,:)],'--r')
            ylabel('other T shock');
            xlabel('Quarters');
            %axis([1 10 -0.4 0.6]);
            subplot(3,5,12)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+pndx,:)],'-k',tt,[nan,-Resplow(2*nvars+pndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+pndx,:)],'--r')
            xlabel('Quarters');
            subplot(3,5,13)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+hndx,:)],'-k',tt,[nan,-Resplow(2*nvars+hndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+hndx,:)],'--r')
            xlabel('Quarters');
            %axis([1 10 -1 0.5]);
            subplot(3,5,14);
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+relendx,:)],'-k',tt,[nan,-Resplow(2*nvars+relendx,:)],'--r',tt,[nan,-Resphigh(2*nvars+relendx,:)],'--r')
            subplot(3,5,15);
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+relsndx,:)],'-k',tt,[nan,-Resplow(2*nvars+relsndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+relsndx,:)],'--r')
            xlabel('Quarters')
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 8])
            print('-depsc2',['SBT_tri_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            
        elseif spec == 28
            figure(1);
            subplot(3,5,1)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
            axis([1 10 -0.5 0.5]);
            title('premium');
            ylabel('SBT shock');
            subplot(3,5,2); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+relendx,:)],'-k',tt,[nan,Resplow(nvars+relendx,:)],'--r',tt,[nan,Resphigh(nvars+relendx,:)],'--r')
            title('rel. hours')
            axis([1 10 -1.2 1.2]);
            subplot(3,5,3)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pndx,:)],'-k',tt,[nan,Resplow(nvars+pndx,:)],'--r',tt,[nan,Resphigh(nvars+pndx,:)],'--r')
            title('prod.');
            axis([1 10 -0.5 0.6]);
            subplot(3,5,4)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hndx,:)],'-k',tt,[nan,Resplow(nvars+hndx,:)],'--r',tt,[nan,Resphigh(nvars+hndx,:)],'--r')
            title('tot. hours');
            axis([1 10 -1 0.5]);
            subplot(3,5,5)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(nvars+pricendx,:)],'-k',tt,[nan,-Resplow(nvars+pricendx,:)],'--r',tt,[nan,-Resphigh(nvars+pricendx,:)],'--r')
            title('price');
            axis([1 10 -1 0.7]);
            subplot(3,5,6)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+spndx,:)],'-k',tt,[nan,-Resplow(2*nvars+spndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+spndx,:)],'--r')
            ylabel('other T shock');
            axis([1 10 -0.5 0.5]);
            subplot(3,5,7);
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+relendx,:)],'-k',tt,[nan,-Resplow(2*nvars+relendx,:)],'--r',tt,[nan,-Resphigh(2*nvars+relendx,:)],'--r')
            axis([1 10 -1.2 1.2]);
            subplot(3,5,8)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+pndx,:)],'-k',tt,[nan,-Resplow(2*nvars+pndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+pndx,:)],'--r')
            axis([1 10 -0.5 0.6]);
            subplot(3,5,9)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+hndx,:)],'-k',tt,[nan,-Resplow(2*nvars+hndx,:)],'--r',tt,[nan,-Resphigh(2*nvars+hndx,:)],'--r')
            axis([1 10 -1 0.5]);
            subplot(3,5,10)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(2*nvars+pricendx,:)],'-k',tt,[nan,Resplow(2*nvars+pricendx,:)],'--r',tt,[nan,Resphigh(2*nvars+pricendx,:)],'--r')
            axis([1 10 -1 0.7]);
            subplot(3,5,11)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(spndx,:)],'-k',tt,[nan,-Resplow(spndx,:)],'--r',tt,[nan,-Resphigh(spndx,:)],'--r')
            ylabel('supply shock');
            xlabel('Quarters');
            axis([1 10 -0.5 0.5]);
            subplot(3,5,12)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(relendx,:)],'-k',tt,[nan,-Resplow(relendx,:)],'--r',tt,[nan,Resphigh(relendx,:)],'--r')
            axis([1 10 -1.2 1.2]);
            xlabel('Quarters');
            subplot(3,5,13);
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(pndx,:)],'-k',tt,[nan,-Resplow(pndx,:)],'--r',tt,[nan,-Resphigh(pndx,:)],'--r')
            axis([1 10 -0.5 0.6]);
            xlabel('Quarters');
            subplot(3,5,14)
            plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(hndx,:)],'-k',tt,[nan,-Resplow(hndx,:)],'--r',tt,[nan,-Resphigh(hndx,:)],'--r')
            xlabel('Quarters');
            axis([1 10 -1 0.5]);
            subplot(3,5,15)
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(pricendx,:)],'-k',tt,[nan,Resplow(pricendx,:)],'--r',tt,[nan,Resphigh(pricendx,:)],'--r')
            axis([1 10 -1 0.7]);
            xlabel('Quarters');
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 20 12])
            if sample == 3 || sample == 4
                display('not saved: choose other sample');
            else
                print('-depsc2',['SBT_triprice_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
            end
        end
    else % do_sign == 0
        if spec == 6 % for Figure 1!!!
            figure(1);
            subplot(1,3,1);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(spndx,:)],'-k',tt,[nan,Resplow(spndx,:)],'--r',tt,[nan,Resphigh(spndx,:)],'--r')
            ylabel('SBT+UBT shock');
            axis([0 10 -0.6 0.7]);
            subplot(1,3,2);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(pndx,:)],'-k',tt,[nan,Resplow(pndx,:)],'--r',tt,[nan,Resphigh(pndx,:)],'--r')
            axis([0 10 -0.6 0.7]);
            subplot(1,3,3);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(hndx,:)],'-k',tt,[nan,Resplow(hndx,:)],'--r',tt,[nan,Resphigh(hndx,:)],'--r')
            axis([0 10 -0.8 0.7]);
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 4])
            print('-depsc2',['SBT_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_SBT.eps'])
            
            figure(2);
            subplot(1,3,1);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+spndx,:)],'-k',tt,[nan,Resplow(nvars+spndx,:)],'--r',tt,[nan,Resphigh(nvars+spndx,:)],'--r')
            ylabel('SNT shock');
            xlabel('Quarters');
            axis([0 10 -0.6 0.7]);
            subplot(1,3,2);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pndx,:)],'-k',tt,[nan,Resplow(nvars+pndx,:)],'--r',tt,[nan,Resphigh(nvars+pndx,:)],'--r')
            xlabel('Quarters');
            axis([0 10 -0.6 0.7]);
            subplot(1,3,3);
            plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+hndx,:)],'-k',tt,[nan,Resplow(nvars+hndx,:)],'--r',tt,[nan,Resphigh(nvars+hndx,:)],'--r')
            xlabel('Quarters');
            axis([0 10 -0.8 0.7]);
            set(gcf,'paperunits','centimeters','PaperPosition',[0 0 13 4])
            print('-depsc2',['SBT_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_SNT.eps'])
        else
            display('no graphs for this option')
        end
    end
    
    % joint identification
elseif spec == 32
    
    if do_sign == 2
        
        figure(1);
        tt = (0:HORIZON_impresp);
        zerobar = zeros(1,HORIZON_impresp);
        subplot(4,5,1); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(pricendx,:)],'-k',tt,[nan,-Resplow(pricendx,:)],'--r',tt,[nan,-Resphigh(pricendx,:)],'--r')
        title('price');
        ylabel('IBT shock');
        axis([0 10 -1.6 0.2]);
        subplot(4,5,2); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(spndx,:)],'-k',tt,[nan,Resplow(spndx,:)],'--r',tt,[nan,Resphigh(spndx,:)],'--r')
        title('premium')
        axis([0 10 -0.5 0.5]);
        subplot(4,5,3); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(relendx,:)],'-k',tt,[nan,Resplow(relendx,:)],'--r',tt,[nan,Resphigh(relendx,:)],'--r')
        title('rel. hours')
        axis([0 10 -1.2 1.2]);
        subplot(4,5,4); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(pndx,:)],'-k',tt,[nan,Resplow(pndx,:)],'--r',tt,[nan,Resphigh(pndx,:)],'--r')
        title('prod.')
        axis([0 10 -0.5 0.6]);
        subplot(4,5,5); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(hndx,:)],'-k',tt,[nan,Resplow(hndx,:)],'--r',tt,[nan,Resphigh(hndx,:)],'--r')
        title('tot.hours');
        axis([0 10 -1 0.7]);
        subplot(4,5,6); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(2*nvars+pricendx,:)],'-k',tt,[nan,-Resplow(nvars+pricendx,:)],'--r',tt,[nan,-Resphigh(nvars+pricendx,:)],'--r')
        ylabel('SBT shock');
        axis([0 10 -1.6 0.2]);
        subplot(4,5,7); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(2*nvars+spndx,:)],'-k',tt,[nan,Resplow(2*nvars+spndx,:)],'--r',tt,[nan,Resphigh(2*nvars+spndx,:)],'--r')
        axis([0 10 -0.5 0.5]);
        subplot(4,5,8); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(2*nvars+relendx,:)],'-k',tt,[nan,Resplow(2*nvars+relendx,:)],'--r',tt,[nan,Resphigh(2*nvars+relendx,:)],'--r')
        axis([0 10 -1.2 1.2]);
        subplot(4,5,9); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(2*nvars+pndx,:)],'-k',tt,[nan,Resplow(2*nvars+pndx,:)],'--r',tt,[nan,Resphigh(2*nvars+pndx,:)],'--r')
        axis([0 10 -0.5 0.6]);
        subplot(4,5,10); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(2*nvars+hndx,:)],'-k',tt,[nan,Resplow(2*nvars+hndx,:)],'--r',tt,[nan,Resphigh(2*nvars+hndx,:)],'--r')
        axis([0 10 -1 0.7]);
        subplot(4,5,11); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(3*nvars+pricendx,:)],'-k',tt,[nan,Resplow(2*nvars+pricendx,:)],'--r',tt,[nan,Resphigh(2*nvars+pricendx,:)],'--r')
        ylabel('other T shock');
        axis([0 10 -1.6 0.2]);
        subplot(4,5,12); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(3*nvars+spndx,:)],'-k',tt,[nan,-Resplow(3*nvars+spndx,:)],'--r',tt,[nan,-Resphigh(3*nvars+spndx,:)],'--r')
        axis([0 10 -0.5 0.5]);
        subplot(4,5,13); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(3*nvars+relendx,:)],'-k',tt,[nan,-Resplow(3*nvars+relendx,:)],'--r',tt,[nan,-Resphigh(3*nvars+relendx,:)],'--r')
        axis([0 10 -1.2 1.2]);
        subplot(4,5,14); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(3*nvars+pndx,:)],'-k',tt,[nan,-Resplow(3*nvars+pndx,:)],'--r',tt,[nan,-Resphigh(3*nvars+pndx,:)],'--r')
        axis([0 10 -0.5 0.6]);
        subplot(4,5,15); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(3*nvars+hndx,:)],'-k',tt,[nan,-Resplow(3*nvars+hndx,:)],'--r',tt,[nan,-Resphigh(3*nvars+hndx,:)],'--r')
        axis([0 10 -1 0.7]);
        subplot(4,5,16); plot(tt,[0,zerobar],'-k',tt,[nan,Respmed(nvars+pricendx,:)],'-k',tt,[nan,Resplow(nvars+pricendx,:)],'--r',tt,[nan,Resphigh(nvars+pricendx,:)],'--r')
        ylabel('supply shock');
        xlabel('Quarters');
        axis([0 10 -1.6 0.2]);
        subplot(4,5,17); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(nvars+spndx,:)],'-k',tt,[nan,-Resplow(nvars+spndx,:)],'--r',tt,[nan,-Resphigh(nvars+spndx,:)],'--r')
        axis([0 10 -0.5 0.5]);
        xlabel('Quarters');
        subplot(4,5,18); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(nvars+relendx,:)],'-k',tt,[nan,-Resplow(nvars+relendx,:)],'--r',tt,[nan,-Resphigh(nvars+relendx,:)],'--r')
        axis([0 10 -1.2 1.2]);
        xlabel('Quarters');
        subplot(4,5,19); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(nvars+pndx,:)],'-k',tt,[nan,-Resplow(nvars+pndx,:)],'--r',tt,[nan,-Resphigh(nvars+pndx,:)],'--r')
        axis([0 10 -0.5 0.6]);
        xlabel('Quarters');
        subplot(4,5,20); plot(tt,[0,zerobar],'-k',tt,[nan,-Respmed(nvars+hndx,:)],'-k',tt,[nan,-Resplow(nvars+hndx,:)],'--r',tt,[nan,-Resphigh(nvars+hndx,:)],'--r')
        axis([0 10 -1 0.7]);
        xlabel('Quarters');
        set(gcf,'paperunits','centimeters','PaperPosition',[0 0 20 12])
        if sample == 3 || sample ==4
            display('not printed')
        else
            print('-depsc2',['jointI_tri_spec_' int2str(spec) '_sample' int2str(sample) '_prem_' int2str(prem_choice) '_prod' int2str(prod_choice) '_price' int2str(price_choice) '_hours' int2str(hours_choice) '_hlev' int2str(do_level) '_break' int2str(do_break) int2str(break_choice) '_est' int2str(prior) int2str(nlagsvar) int2str(decay) '_sign' '.eps'])
        end
        
        
    end
    
    
end
